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Abstract 

We present a quantization of the hydrodynamic model to describe the excitation of plasmons in a single- 
walled carbon nanotube by a fast point charge moving near its surface at an ai^bitraiy angle of incidence. 
Using a two-dimensional electron gas represented by two interacting fluids, which takes into account the 
different nature of the a and vr electrons, we obtain plasmon energies in near-quantitative agreement with 
experiment. Further, the implemented quantization procedure allows us to study the probability of exciting 
various plasmon modes, as well as the stopping force and energy loss spectra of the incident particle. 
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I. INTRODUCTION 



Plasmon excitations in carbon nanotubes (CNTs) continue to attract attention for a variety of 
applications, e.g., in the context of their optical response llMj]. On the other hand, plasmon ex- 
citations are most effectively probed by fast-moving charged particles, such as in electron energy 
loss spectroscopy (EELS) using a scanning transmission electron microscope (STEM). This tech- 
nique has proven to be a powerful tool for investigating the dynamic response of CNTs i^HM and, 
more recently, graphene |8, ^. Moreover, plasmon excitations were observed to play an impor- 
tant role in electronic excitations of CNTs exposed to ion bombardment llioll . On the theoretical 
side, in addition to ab initio calculations [jg, [Uj], simpler models have been also used to study 



plasmon excitations in carbon nano-structures, such as the empirical dielectric tensor [121 and the 
hydrodynamic model lll3n . 

Since carbon nano-structures present physical realizations of one-atom thick layers of an elec- 
tron gas, it is no surprise that a two-dimensional (2D) version of the hydrodynamic model was used 
early on in studying plasmon excitation in CNTs [Il4l-ll6n. We note that a planar 2D hydrodynamic 



In, 



1811 in 1973. This model has subsequently 



model of the electron gas was pioneered by Fetter 
been used to gain qualitative understanding of plasmon excitations in various quasi-2D electronic 
structures, including semiconductor inversion layers, quantum wells, and thin metallic films. Such 
a model assumes that all electrons belong to a single fluid characterized by two parameters only, 
the equilibrium surface density uq and the effective electron mass m*. Owing to its simplicity and 
versatility in handling difficult geometric constraints, multilayered structures, and the presence of 
dielectric environment, the single-fluid version of the 2D hydrodynamic model has been used in a 



significant number of applications of CNTs and FuUerene molecules [[19- 

However, treating the four valence electrons per atom in carbon nano-structures as belonging 
to the same fluid fails to capture significant differences in their bonding. Three of these valence 
electrons are involved in the strong a bonds characterizing the sp"^ hybridization within a layer of 
carbon atoms. At the same time, one electron occupies a weakly bound vr orbital that is largely 
responsible for the conductivity properties of carbon nano-structures. While the vr electron bands 
make graphene a zero-gap semiconductor, the so-called graphene nano-ribbons (GNRs), as well 
as CNTs may exhibit metallic or semiconducting character depending on the symmetry of their 
underlying atomic structure. On the other hand, the a electron bands exhibit a large gap that gives 
rise to a characteristic high-frequency feature in the absorption spectra in various allotropic forms 
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of carbon. 

Consequently, inspired by the empirical model of Cazaux 112 811 . Barton and Eberlein n2% pro- 
posed a two-fluid version of the 2D hydrodynamic model. Their model treats the a and it elec- 
trons as two classical fluids with equilibrium densities = 3 Uat and = Uat (where Uat ~ 38 
nm~^ is graphene's atomic density), respectively, which they superimposed on the surface of a Ceo 
molecule. By introducing an empirical restoring frequency for the a electron fluid, the two-fluid 
2D hydrodynamic model of Barton and Eberlein was shown to give rise to two groups of plas- 
mons. This explained well the two dominant absorption features in the frequency ranges of 5-10 



M and single-wall CNTs (SWCNTs) UM- 



eV and 15-30 eV, in both FuUerene molecules 11291 

On the other hand, Mowbray et al. introduced the quantum, or Fermi pressure in the a and 
TT electron fluids, and reformulated the two-fluid 2D hydrodynamic model without invoking any 
empirical parameters ll32ll . They found that the electrostatic interaction between the fluids gives 
rise to two groups of high- and low-frequency plasmons for both SWCNTs [|32n and multi-wall 
CNTs (MWCNTs) [33] in similar ranges as those obtained by Barton and Eberlein It^. Versatility 
of this version of the two-fluid 2D hydrod yna mic model was further exploited in studying the 
interactions of CNTs with dielectric media [24, 251' as well as in applications to the propagation 
of electromagnetic radiation through CNTs [l36l - l38ll . 

In the present paper, we perform a second quantization of the two-fluid 2D hydrodynamic 
model for electron gas confined on the surface of an infinitely long cylinder representing a free 
SWCNT, and use it to analyze the spectra of plasmon excitations induced by a classical non- 
relativistic charged particle, passing by the nanotube under an arbitrary angle of incidence. This 
approach enables us to study the effects of multiple plasmon excitations by external particle beams, 
similar to the EELS studies of thin metallic films on solid surfaces [39]. Such a study can help 
elucidate which plasmon modes are most effectively excited depending on the incident projectile 
trajectory. Besides obvious applications in the EELS studies of CNTs, the results of this work may 



help to understand the ro 



during ion irradiation [|lOtl 



e of electronic channels in energy deposition in carbon nano-structures 



After describing the model and outlining its second quantization in section [HI we shall discuss 
our results for the plasmon excitation spectra, total energy loss, and the stopping force on the 
external charge in section [HIl Our concluding remarks will be followed by a description of the 
statistical properties of plasmon excitations in appendix[Al a discussion of the relation to the semi- 
classical model in appendix [Bj and an outline of how the results for planar 2D electron gas may be 
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obtained from the formalism for cylindrical geometry in appendix O We shall use the Gaussian 
electrostatic units unless otherwise indicated. 



II. BASIC THEORY 

A. Hydrodynamical two-fluid Hamiltonian 

The total Hamiltonian, for the a and vr electron fluids can be written as 



I dhs n,(r„ t) I ^M^s, t) - 



32 , , 



1/2- 



+ Y 5^ E jjd'r.dh'^'^^^^^^ (1) 



where = {r = R] z} is a position on the surface of a single-wall carbon nanotube of radius R 
aligned with the z-axis of a cylindrical coordinate system with coordinates r = {r, z}, whereas 
Uy, Tpy, , diWA Ky , rc spcc tlvcly, the number density per unit area, the fluid momentum 

field, displacement field, effective mass, and the restoring-force constant for the z/th electron fluid, 
with u = Ti.a. Note that the first term in Eq. ([U) may be interpreted as the kinetic energy of a clas- 
sical fluid V moving at the velocity vl^{ys, t) = Pui^Cs, t)/ml = x,^(rs, t), where the dot stands for 
time derivative. The second term represents the restoring (R) effects on the fluid displacement due 
to electron binding in harmonic approximation, giving rise to restoring frequency uj^r = a/ Ky/ml. 
Further, the potential energy per electron, V = Vg^ + Kxt^ consists of the ground state energy due 
to the positive-ion background, Vgr{rs), and the energy due to the time-dependent external perturb- 
ing potential, Vext(rs7 1)- The next two terms represent, respectively, the Thomas-Fermi (TF) and 
the Dirac's (D) exchange interactions in the local-density approximation for the electron fluid u, 
whereas the last term in Eq. ([T]) represents Coulomb interactions among electrons. 

A few additional comments on the Hamiltonian in Eq. ([T]) may be in order. While restoring 
frequencies ujyj. play the role of empirical parameters that are usually determined by fitting to 
experimental data their physical motivation is to take into account the effects of 

the electronic band structure of carbon nano-structures in a manner analogous to that invoked in 

o 

devising the Drude-Lorentz dielectric function for carbon materials [|40|] . So, invoking the high 
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mobility of the vr electrons in metallic nanotubes and considering that the a electron bands exhibit 
a gap of about 12 eV, we adopt the resto ring frequencies of uj-^r = eV and oj^r = 16 eV, which 
have been used by previous authors [|29l-l3 IN. Next, we note that the TF and D terms are written 
in the Hamiltonian in Eq. ([T]) in a form commensurate with a planar 2D electron gas. This may 
be justified by considering a quasi-free 2D electron gas with the number density per unit area n, 
occupying the surface of a cylinder with radius R, where curvature effects were found negligible 



when Rkp ^ 1, where kp = \phvn is the Fermi wavenumber of the corresponding planar 2D 
electron gas [1341. kill. Finally, we note that a Hamiltonian, which is similar in form to the one in 
Eq. (dl), was used by van Zyl and Zaremba 1421] who included the so-called von Weizsacker, or 
gradient-correction term of the form 



vW 



5m: 



|Vri.(r„t)|p 



(2) 



This provides an approximate correction for the non-local effects in the TF interaction in a planar 
2D electron gas. However, there is considerable uncertainty regarding the exact form of such a 
term. In fact, it was confirmed recently that the von Weizsacker factor A^w in Eq. Q should vanish 
in a strict planar case ll43n . Nevertheless, it has become quite common to use the von Weizsacker 
correction with Avw = 1 in many recent studies of CNTs [2^ ZllSI^], even though the effects 



of curvature on this term are not known at present. Fortunately, the von Weizsacker correction 
was found to only affect plasmon dispersion at very short wavelengths ||34|]. Such wavelengths are 
comparable to the inter-atomic spacing where the hydrodynamic model is likely to break down 
anyway, so we shall neglect such a contribution to the Hamiltonian in Eq. ©. 

We further assume that the interaction with the external particle can be considered as a small 
perturbation. Hence, we can expand the Hamiltonian in Eq. © to the second order with respect 
to the perturbation by declaring V^xt = AKxt and writing \\.v{vs,t) = X5:k^{rs,t) + O(A^) and 
n^{rs,t) = n° + \5ny{Ts,t) + C(A^), where n° is the ground-state electron density in the z/th 
fluid. We can further express the perturbed electron density in terms of the displacement field in 
a manner that will automatically satisfy the linearized continuity equation by writing 5n^{Ts, t) = 
— ?7,°Vj ■ 5x^{rs,t). Here the gradient Vg differentiates only in the directions tangential to the 
surface of the cylinder r = R. Moreover, by restricting consideration to electrostatic phenomena, 
the velocity field will be irrotational. This will be ensured by defining a potential function for the 
displacement field, ^u{^s,t), such that 5xi^(rs,t) = —Vs^u{'^s,t)- We can now write the second- 
order Hamiltonian in terms of the function ^^(rs, t) alone. 



2 r r ^0 ^0 
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Before proceeding, it is worth mentioning that, while quantum treatment generally requires a 



second-order Hamiltonian ll44ll . as in Eq. ([3]), such an approximation for the Hamiltonian yields 
semi-classical equations of motion for the electron gas that correspond to the linearized hydrody- 
namic model. To illustrate this point, we note that one can derive from M2 a Lagrange's equation 
of motion for the potential function ^^{vs, t) as follows 

't + ujlriu - slV^i^ = ^ [Kxt(rs, t) - e<l>i„d(r„ t)] , (4) 

where we have defined the speed of propagation of density disturbances in the z/th fluid, s^, by 

Here ^i„i{r,t) is the induced potential in the system due to polarization of the electron fluids, 
satisfying the Poisson equation in three dimensions, 

V2$M(r, t) = 4ne6{r - R) Yl ^^'U^s, t), (6) 

U=7V,(7 

where we have used the fact that the induced density in the z/th fluid is given by 6ni, = n^JW^^y 
to the first order. It is interesting to note that the effect of the Dirac term in Eq. ([5]) is to reduce 
the adiabatic limit of the speed due to the TF term, which can be written as stfu = vpu/ with 
vfu = (/i/m*)-\/27™° being the Fermi speed for the i^-th fluid in equilibrium. It can be shown that. 



in a hydrodynamic model without the Dirac correction 1 1^ 3.|32|], the high-frequency expression 



for the speed due to the TF interaction should be corrected so that stfu = {'\/^/'^)vfu [I45ll . 



However, including this effect in the present model would require using a relaxation approximation 



for the kinetics of the electron gas, rendering second quantization of the model unsuitable [|45|]. 

We make further progress by seeking a Fourier series representation of the potential function 
^ui'i^s, t) = ^ui'^, z, t) for its dependence on coordinates on the nanotube's surface. This will give 
rise to a change in variables {(/?, z} 1— {m, k} for all quantities of interest, with m and k defining 
the modes of plasmon propagation around the nanotube's circumference and in its axial direction. 



respectively. We initially assume that the nanotube has a finite length L, so that its total area is 
A = 2'kRL and A; is a discrete variable, taking values k = with i being an integer. Thus, 
defining the Fourier coefficients as 

-L/2 



iv{m, k,t) = R dip 

J —IT 

we may write 

^ oo 

Urs,t) = - J2 E 



L/2 



J,rrnp+ikz 



(7) 



dk 
— < 
27r 



^irmp+ikz 



Um,k,t),(^) 



m.=— oo k 

where the rightmost expression involves a Fourier transform in the axial direction that arises in the 
limit L — 7- oo, when k becomes a continuous variable. We note that, while the present formalism is 
most clearly developed by using summation over k for a nanotube of finite length, our final results 
will be computed by integration over k, corresponding to the limit of an infinitely long nanotube. 
In this way we neglect any end effects on plasmon excitation spectra, which may be justified by 
invoking a typically high aspect ratio for CNTs, L/R ^ 1. 
One can now express M2 in Eq. © as 



m,k ^ 



"'mk ^ 



m.k 



(9) 



with the unperturbed part of the Hamiltonian given by 

.2 



h 



(0) 
mk 



E 



n 



0^ 
" 2 



m 



i?2 



(10) 



2 / 2 \ 

+ ^Rgmk f ^ + j I] nlnl,Ct' 
and the interacting part by 



(11) 



where we have suppressed the dependencies on (m, k, t) in and its complex conjugate Here, 
Qmk = ^iTrlmilklR) Km{\k\R) is the Green's function for the Poisson equation in cylindrical coor- 
dinates evaluated at the nanotube's surface, with and Km being the modified Bessel functions of 
integer order m, of the first and second kind, respectively. Moreover, Vcxt = Kxt(-R; fn, k, t) is the 
Fourier transform with respect to coordinates of the external potential Kxtlfs, t) = Kxt(-R; z, t), 
evaluated on the nanotube's surface. In particular, for a point charge eZ moving on a trajectory 
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ro(t) = {ro{t),(po{t), zo{t)} that remains external to the nanotube at all times, ro(t) > R, we 
obtain 

V^^R; m, k, t) = -e'Z Rg^^iR, ro(t)) e-*-'^°W-^'=^«W, (12) 

where g.^k{R,ro) = 47r/„(|A;|i?)/C;„(|/i;|ro). 

We now proceed to diagonalize /i^^ in the subspace of interactinga and vr fluids for fixed m 



and k. (The procedure is similar to that followed by Gorokhov et al. BOh in their treatment of C, 



-60 



plasmon excitation with a quantized two-fluid model). To this end, we define the factors 



— for a electrons 



/. = = < (13) 

7— for vr electrons, 

where uq = + n1 is the total areal density of valence electrons in graphene, with relative 



weights 3/4 and 1/4 corresponding to the sp^ hybridization [|40n . and is a suitably defined mean 



effective electron mass (see below). Thus, we can rewrite Eq. (flOl ) as 



h^'l = [/<. {I'w + u^lu') + (le.p + c^^ie.p) + ^^^^^ (a. + ce. 



(14) 



where we have defined plasma frequencies of non-interacting fluids, cOa- and cj^, by 

m* J \ R^ 



= ^1 +{sl + e'Rg^S\ i'^ + k']. (15) 



while their coupling is defined by 



^' = e^Rgm,\l^^['^,+k^]. (16) 

One can diagonalize the Hamiltonian by substituting 

V^To- ^o- = ^1 cos a — /I2 sin a, (17) 



A/X-,^7r = ^1 sina + cosa, (18) 

into Eq. (fT4l) and by choosing the angle a so that oscillations with amplitudes Ai and are decou- 
pled. This is achieved when cot (2a) = ^^^f^r^, so that the non-interacting part of the Hamiltonian 
can be written as a sum of decoupled oscillators. 



h 



(0) 
mk 



i=i 



(19) 



where we have restored the dependencies on (m, fc, t) and defined the eigen-frequencies, Ujmk > 0, 
of the decoupled oscillators by 



^ jmk 



± 



+ A4 



(20) 



with j = 1,2 for the positive and negative signs, respectively. 

For the interacting Hamiltonian, we substitute Eqs. (fTTI) and (fTSi) into Eq. (fTT)) and obtain 



h 



(int) 
mk 



1 n" 



2v^ 

1 



2v^ 



{Ai cos a — ^2 sin a) V*^^ + c.a 
{Ai sin a + A2 cos a) V^*^ + c.c. 



(21) 



which can be written in a more compact form if we define m^, to be a weighted harmonic mean of 
the effective masses in the a and vr fluids, 



as follows 

where we have restored the dependencies on (m, /c, t). Here 

-cos(a-/3), j = 1, 
sin(a-/3), j = 2, 

where the angles a and (3 are defined, respectively, by 



(22) 



(23) 



jmk 



(24) 



cos a 



sm a 



V2' 



1 ± 



UJ 



Imk 



(25) 



and 



cos f3 



sin /3 



2mk 



(26) 



It is worth mentioning that a classical, single-fluid model is recovered by letting both u^r — )■ 
and s,^ — > in Eq. ( fTSl ). so that Eq. (|20l ) then gives a = /3 with u;2mfc = and Wi^fe corresponding 
to plasma frequency of a 2D electron gas with surface density no = n° + n° and an effective mass 



defined via Eq.daa) 12 



3, 
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B. Quantization 



Once we have the decoupled Hamiltonian to describe the system, it is useful to apply a quan- 
tum treatment in terms of creation and annihilation operators. This provides a simple and clear 
way to describe the excitation of plasmons (oscillators) due to the interaction with the external 
particle. Following the quantization procedure described by Arista and Fuentes n46\ . we assign to 
the coefficients Ajm^ creation and annihilation operators as follows 



jmk '- 



O-jruki^) + ajmk{t) 



(27) 



where the coefficient ■jjmk is to be determined. The operators ajmk and dj^^, satisfy the usual 
commutation relations 

t 



ttj'm'k', Ctjmk 



^jj' ^kk' ^mm' i 

while their time dependence, djmk{t) = a jmk exp {-iujmkt) , gives ajmk{t) = ^^jmk(^jmk 

(t), SO 

that 



Ajmk (^) 



• '~1jmk 



^^]mkif) - ajmkit) 



(28) 



With these relations, the (quantized) non-interacting Hamiltonian of the mode (m, k) can be 
written as 



h 



(0) 
mk 



m 



no- 



2 / -t 
2 Ijmk \ '^jmk^jmk 



(29) 



giving for the full non-interacting Hamiltonian (see Eq. (9)) 



m,k 

= ^ hUjmk {a'jmk 



(0) 
mk 



(30) 



This procedure gives the coefficient 7jmfc to be used in Eq. (1271) as 

''ijmk ~ 



mk 



For the (quantized) interacting Hamiltonian of the mode (m, k) we obtain 



(31) 



h 



(int) _ nQ 



mk 



''Ijmk 



2 ^UJjjjik 



D 



jmk 



O'lmki't) + ajmk{t) K*xt(^; t) + h. 



(32) 
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so that the total interacting Hamiltonian reads 



(33) 



with 



Vcxtl-R; k,t) 



(34) 



We note that the interacting Hamiltonian in Eq. (l33l) is given by a superposition of the linear 
displacements of plasmons as quantum oscillators, which is a consequence of using the second- 
order approximation for the full Hamiltonian, given in Eq. (|3]) Hm . 

We further follow the formalism presented by Arista and Fuentes 1461 and obtain the average 
number of excited plasmons in a given mode (j, m, k) as 

1 



to 



where to ^ ^ oo and t — )■ +oo. Inserting Eq. (|34l) in Eq. (I35|) . we obtain 



jmk 



Kxt fc, 



jmk ) 



where 



jmk 



r)2 mi A_ U2 
^jmk r2 ^ 



(35) 



(36) 



(37) 



A 2hm^ oOjmk 

and Kxt(-R; k, u) is the Fourier transform with respect to time of the function V^xt(-R; 'm, k, t), 

/"OO 

VUR\m.k,uj)= dte'^XAR:m,k,t). (38) 



The time integral in Eq. (1381) can be solved analytically by using Eq. (112]) for the case of a 
straight-line trajectory with constant velocity, given by ro(t) = {ro{t), ipo{t), ZQ{t)} with 



roit) 


= \/rLn + vlt^^ 




= arctan f^'-M 


\ '"min / 


Zoit) 


= V\\t, 



where rmin > -R is the shortest distance between the trajectory and the nanotube axis, and v_i 
and are the components of the incident particle's velocity in the directions perpendicular and 
parallel to the nanotube's axis, respectively. Thus, one obtains n3m 



VeniR;m,k,u) = -iTie^ Z RImi\k\R)}Cmik, uj) , 



(39) 
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where 



) 



m 



(40) 



= vr 



^(a;-^l)2 + (A;i;x)2 



Here, the velocity components may be expressed in terms of the total projectile speed, v, and 
the incident angle, Q, relative to the nanotube axis as = t>sin6' and v\\ = vcosO. Without 
any loss of generality, one may adopt the range of angles < 9 < 7r/2, so that the direction of 
projectile motion is defined by positive values v± > and v^^ > 0, corresponding to the directions 
of increasing variables ip and z, respectively. As a consequence, one can then infer from Eq. (|40l) 
that the propagation direction of plasmon modes (m, k), which is commensurate with the direction 
of projectile motion, will be given by positive values, m > and A; > 0. 

Finally, using the fact that the plasma excitations by an external particle are represented by 
independent quantum oscillators of the mode (j, m, /c), and that the probability distribution of 
exciting such an oscillator to the A^th state is Poissonian defined by the mean fj.jmk, one can 
deduce several statistical properties of practical interest. For example, it is shown in appendix lAl 
that the probability distribution of plasma excitations in the mode with fixed j and m for arbitrary 
k is also Poissonian, defined by the mean 



Moreover, we shall derive in appendix |A] an expression for the probability density for energy loss, 
P{e) in Eq. (IA7I) . as well as the total energy loss, AEioss in Eq. (IA8I) . which can be compared with 
experiments or other theoretical approaches. 

A few comments may be in order about the projectile trajectory in relation to plasmon excita- 
tions of the nanotube. The straight-line trajectory presents a good approximation for a projectile 
with sufficiently high momentum, so that the effects of the recoil and trajectory bending may be 
neglected. On the other hand, a projectile moving at very high speeds would require inclusion 
of the retardation effects in the electron response of the CNT, which is beyond the scope of the 
present work. Moreover, our neglect of the end effects on plasmon excitation may be justified 
when the projectile moves on non-parallel trajectories with angles tan 9 ^ R/ L. In practice, this 
condition is satisfied for angles 9 > 1° owing to the high aspect ratio of CNTs. 

The case of parallel trajectory with 9 = deserves special attention. In that case, one has 
to assume that the nanotube length L is finite, but still large enough to allow the neglect of end 
effects on plasmon excitation, as indicated in the limiting form of Eq. ([8]). On the other hand. 




(41) 



12 



it is rigorously shown in the last paragraph of appendix |B] that, when v± = 0, the number of 
excited plasmons in the mode (j, m, k) becomes proportional to L6(ujmk — ki'n), and Eq. (|36l) 
must be then evaluated by using Eq. (IB 161 ). rather than Eqs. (|39l ) and (|40l ). If L may also be 
considered short enough, one can further define the stopping force, or the energy loss per unit path 
length, S = AEioss/ L, giving Eq. (IB 171 ) for a projectile moving on a parallel trajectory outside the 
nanotube. That result is extended in a straightforward manner in Eq. (IB 181) for a particle channeled 



inside the nanotube 124 



2511 ■ In either case, the quantity S may be considered independent of L 
if the total energy loss A^Eioss in a finite-length nanotube is much smaller than the initial kinetic 
energy of the particle on a parallel trajectory. 

III. RESULTS AND DISCUSSION 

We first compare in Fig. HJa) the effects of restoring (R), Thomas-Fermi (TF) and Dirac (D) 
interactions in Eq. ([!]) on the dispersion relations ujjm{k), defined in Eq. (|20|) . for modes j = 1,2 
with m = in a CNT with radius i? = 7 A, assuming that the effective masses are all equal to 
the free electron mass. We note that setting all these interactions to zero causes the two fluids 
to collapse into a single classical electron fluid (SF) with surface density uq = Ans_t, giving rise 
to the cr + 7r electron plasma oscillations [I2O]. Introducing either the R interaction or the TF 
(with or without the D) interaction terms causes splitting of the plasmon dispersions into the high- 
frequency ij = 1 or a + n) branch, which is very close to the single-fluid branch, and the low 
frequency (j = 2 or tt) branch. The main difference between the splitting of plasmon branches due 
to the R interaction as opposed to the TF(D) interaction is that, in the latter case, both the upper 
and the lower plasmon dispersions tUjok vanish as — > 0, while the tt branch is markedly acoustic. 
On the other hand, inclusion of the R interaction renders the tt plasmon branch quasi-acoustic, 
with an almost constant slope for A; > 0.1 A~^, in close agreement with the experimental data 
[4,^7]. It is also seen in Fig. [TJa) that the inclusion of the D interaction both lowers the upper 
branch to some extent and reduces the slope of the lower branch, as expected from Eq. (|5]). 

Since in the rest of this work we use a "complete" two-fluid hydrodynamic model, which 
includes restoring (R), Thomas-Fermi (TF) and Dirac (D) interactions, we compare in Fig. [TJb) 
its modes with m = 0, 1, 2, 3, and 4 for both the upper (7 = 1 or cr + tt) and lower (j = 2 



or tt) plasmon branches with the experimental data [|4 , ^ |7|] . This semi-quantitative agreement 
with experiment, particularly for the narrow group of lower plasmon branches, indicates that the 
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MA] 



FIG. 1: Plasmon energies huo in eV versus longitudinal momentum in A ^ of a SWCNT with radius 
R = 7 A for (a) angular momentum m = from single-fluid (SF) 12011. two-fluid Thomas-Fermi (TF) 



1128 



350, two-fluid Thomas-Fermi-Dirac (TFD), restoring frequency (R) ll29ll . two-fluid Thomas-Fermi with 



restoring frequency (TFR), and two fluid Thomas-Fermi-Dirac with restoring frequency (TFDR) models, 
and (b) two fluid Thomas-Fermi-Dirac with restoring frequency model for angular momenta m = 0, 1, 2, 3, 
and 4. Experimental results from Ref. M\ (circles), Ref. (diamonds), and Ref. 171] (squares) are provided 
for comparison. 



inclusion of both the R and D interactions improves the present, complete hydrodynamic model in 
comparison to the previous two-fluid model, which only included the TF interaction [|32ll . 

In the remaining calculations we consider a particle with Z = ±1 (proton or electron), pass- 
ing by a SWCNT of radius i? = 7 A at a minimum distance rmin = 10.5 A (unless indicated 
otherwise). 

Figures [2] and |3] show the average number of excited plasmons, fXjmk, as a function of k for 
|m| = 0, 1, 2, 3, and 4 with j = I (a,c) and j = 2 (b,d), for two total speeds, v = 5 a.u. (a,b) 
and V = 10 a.u. (c,d), and for two incident angles, 6 = 90° (Fig. |2]) and 9 = 45° (Fig. |3]). 
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FIG. 2: Average number of excited plasmons 2^^^jmk as a function of longitudinal momentum k in 

on a SWCNT with R = 7 A for (a,c) o" + vr modes (j = 1) and (b,d) tt modes (j = 2) with angular momenta 

|m| = ( ), 1 ( ), 2 ( ), 3 (- • -), and 4 (- •• -), shown by black lines when m > and orange 

lines when m < 0, due to a proton with rmin = 10.5 A, incident at angle relative to the SWCNT axis of 
9 = 90° with speed (a,b) v = 5 a.u. and (c,d) f = 10 a.u. 

Note that, for normal incidence, /i^mfc is an even function of wavenumber, so only positive k 
values are shown in Fig. [2l We observe that, for low \k\ values, the low-frequency branches 
(j = 2) dominate over the high-energy modes (j = 1) in both Fig. [2] and Fig. [3l but this trend 
seems to be reversed for increasing values of increasing speed, and oblique incidence. In 
particular, the singular behavior of the mode j = 2, m = as — t- in Figs. [2] and [3] stems 
from the way how its dispersion relation vanishes in the long wavelength limit in Fig. [H which 
may be approximated by a one-dimensional plasmon dispersion for the vr electron fluid, uj20k ~ 
k^JA^xRe'^{nl/ml) ln(1.123/fci?) for A; < 0.1 A"^ [41|]. At oblique incidence the iijrnk curves 
are no longer even functions of k, and the observed dominance in Fig.[3]of the values with A; > 
indicates that the preferential direction for propagation of the plasmon coincides with the direction 
of the longitudinal velocity component of the projectile, v\\ > 0. In addition, one notices in Fig. [2] 
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FIG. 3: Average number of excited plasmons 2^^^jmk as a function of longitudinal momentum k in 

on a SWCNT with R = 7 A for (a,c) o" + vr modes (j = 1) and (b,d) tt modes (j = 2) with angular momenta 

|m| = ( ), 1 ( ), 2 ( ), 3 (- • -), and 4 (- •• -), shown by black lines when m > and orange 

lines when m < 0, due to a proton with rmin = 10.5 A, incident at angle relative to the SWCNT axis of 
9 = 45° with speed (a,b) v = 5 a.u. and (c,d) f = 10 a.u. 

that the modes with m > (black lines) are generally excited with much higher probabilities than 
the modes with m < (orange lines) indicating that, for normal incidence, plasmon waves tend to 
propagate around the nanotube circumference in the direction of projectile motion. However, the 
situation is reversed for oblique incidence, so that the modes with m < tend to dominate over 
those with m > for increasing k > values in Fig. |3l 

Figure |4] shows the dependence of the average number of excited plasmons, /ij^, which is 
obtained after integration of the curves in Figs. [2] and [3] over k, on the total particle speed v for 
9 = 90° (a,b) and 6 = 45° (c,d), for the branches j = 1 (a,c) and j = 2 (b,d). One can say that, 
globally, the mode j = 2, m = dominates for any speed or direction of the incident particle, 
displayed in Fig. HI but we also observe that the relative contributions of the other modes change 
with increasing v. For higher speeds of the incident particle, the relative importance of the modes 
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FIG. 4: Average number of excited plasmons ^jm, defined in Eq. ((411) . on a SWCNT witli = 7 A for 

(a,c) (T + TT modes (j = 1) and tt modes {j = 2) with angular momenta |?n| = ( ), 1 ( ), 2 ( — 

-), 3 (- • -), and 4 (- •• -), siiown by black lines when in > and orange lines when m < as a function 
of the total speed v in a.u. of a proton with rmin = 10.5 A at an angle relative to the SWCNT axis of (a,b) 
e = 90° and (c,d) = 45°. 

with j = 1 increases, especially for m = and m. = 1. Moreover, one notices that the modes 
with m < are generally suppressed compared to the modes with m > at lower speeds, but they 
tend to regain some importance as the speed increases. We also notice that there exist velocity 
thresholds for the excitation of the different non-acoustic modes (all j = I modes, and the j = 2 
modes with m ^ 0). We remark that finite values of the average number fijm for the mode j = 2, 
m = 0, which extend down to low speeds in Fig. HI are likely a consequence of the singular 
behavior of the corresponding number /ijmfc as A; — )■ seen in Figs. [2]and|3l but they are not in 
violation of the conservation of energy for slowly moving projectiles, as documented by the total 



energy loss at low speeds in Fig.l7] [|48ll 



The dependence of fXjm on the incidence angle 9 is displayed in Fig.[5]for the modes with j = 1 
(a,c) and j = 2 (b,d) for two total speeds: v = 5 a.u. (a,b) and v = 10 a.u. (c,d). In general, the 
relative importance of various modes is similar to that seen in Fig. |4] exhibiting dominance of the 
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FIG. 5: Average number of excited plasmons fijm, defined in Eq. ((411) . on a SWCNT witli = 7 A for 

(a,c) (T + TT modes (j = 1) and vr modes (j = 2) with angular momenta |?n| = ( ), 1 ( ), 2 ( — 

-), 3 (- • -), and 4 (- •• -), shown by black lines when ?n > and orange lines when m < as a function 
of the angle relative to the SWCNT axis 6 in degrees for a proton with rmin = 10.5 A at the total speed of 
(a,b) V = 5 a.u. and (c,d) v = \0 a.u. 

mode j = 2, m = 0. One notices that the modes with m < 0, which are suppressed at finite angles 
of incidence, become equal to the modes with m > in the limit of grazing incidence, 6* — )■ 0, 
as expected from symmetry considerations. Furthermore, the overall increase in the number of 
all excited modes with decreasing angle 9 may be ascribed to increasing interaction time with an 
otherwise infinite nanotube, which should become proportional to cosec 9 for angles 9 smaller than 
several degrees. In particular, it is interesting to note that, while in Figs. [21 [3] and |4] the average 
numbers of excited plasmon modes are generally small, and even the most dominant mode j = 2, 
m = hardly exceeds unity, the number of plasmons excited in the mode j = 2, m = is seen in 
Fig.[5]to greatly exceed unity for, e.g., 9 < 5°. 

From the results of previous section it is possible to obtain a statistical description of energy 
loss of the incident particle (see appendix lAl). One of the most useful quantities is the probability 
density, P{e), for losing a given amount of energy, e, which can be related to the energy loss 
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FIG. 6: Probability density for proton energy loss e in eV, expressed as eP{e), for proton trajectory passing 
near a SWCNT with R = 7 A at the minimum separation of = 10.5 A with an angle relative to the 
SWCNT axis of (a) 6 = 90° and (b) 9 = 45°, at the total speed of v = \, 2, 5, 10, 20, and 50 a.u.. The 
spectra are vertically shifted (in steps of 0.2) for the sake of clarity. The position of the /c — )• plasmon 
modes Ujmo are provided for comparison. 

spectra with peaks corresponding to the excitation of various plasmon modes. In Fig. |6]we show 
the product sP{e) for several incident speeds of a particle travelling perpendicular (a) and with 
an inclination angle of 45° (b) with respect to the tube axis. We can see two distinct groups of 
peaks, indicating excitation of the two branches of plasmons, vr and cr + vr, shown in Fig. d] At low 
speeds (f ~ 1 — 2 a.u.), only the low-energy plasmons are excited [|49|] . which is commensurate 
with the demonstrated ability of the low-energy EELS technique to probe the vr plasmon excitation 
in graphene [|9|]. On the other hand, as the incident speed increases, we notice in Fig. 6 excitation 
of both a + n and vr plasmons, as observed in the high-energy EELS experiments [l6|-l8|]. Even 
though we do not pursue here the problem of plasmon excitations at relativistic projectile speeds, 
one may infer from Fig.[6]that the excitation probabilities of both groups of plasmons decrease as 
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FIG. 7: Total energy loss Ai?ioss in eV versus speed v in a.u. for proton trajectory passing near a SWCNT 
with R = l k with an angle relative to the SWCNT axis of (a) B = 90° and (b) 6 = 45° , at the minimum 
sepai-ation of rmin = V.5 A (top), 8.0 A, 8.5 A, 9.0 A, 9.5 A, 10.0 A, and 10.5 A (bottom). 

the speed exceeds a value on the order of 20 a.u., with the low-energy plasmon peaks being almost 
diminished and the high-energy peaks still visible at v ~ 50 a.u.. 

Looking into details of the spectra in Fig. [6l one notices that, as the speed increases, there are 
changes not only in the overall intensity of the peaks, but also in the relative weight of each mode. 
At lower speeds, modes m = 2, 3, . . . give large contributions in both branches, while at higher 
speeds their weights are suppressed. It is interesting to notice that the mode j = 1, m = only 
appears at the speeds v > 10 a.u. as a broad feature just above the restoring frequency of 16 eV of 
the a fluid. When the incidence is oblique, we observe broadening of all peaks, especially at low 
speeds, which is due to an increase of the high-/;; contributions to the excitation of various plasmon 
modes (c.f. Figs. [2] and (3]) . 

By integrating the curves shown in Fig.[6l one can obtain the total energy loss, AEioss, suffered 
by the incident particle. This is shown in Fig. |7] versus the particle speed v for two angles of 
incidence, (a) 9 = 90° and (b) 9 = 45°, and for several values of rmm- For the smallest rmin, the 
energy loss presents a large maximum around the speed of v = 2.5 a.u., and a smaller peak at a 
speed V < 1 a.u. The physical reason for the appearance of two peaks is likely due to the relatively 
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V [a.u.] 

FIG. 8: Stopping force in eV/A versus speed v in a.u. for a proton travelling parallel to the axis of a SWCNT 
with R = l k{a) inside the SWCNT at rmin = (bottom), 0.5 A, 1.0 A, . . . , 6.0 A (top), and (b) outside 
the SWCNT at r^in = 8.0 A (top), 8.5 A, ... , 10.5 A (bottom). 

broad gap between the low-energy and high-energy groups of peaks seen in Fig.[6l As expected, 
the two peaks in Fig. |7] decrease in magnitude as the minimum distance rmin increases, and they 
broaden and move towards higher speeds, so that the smaller peak turns into a shoulder around 
V \ — 2 a.u.. It is also interesting to notice that there is a threshold for energy loss on the order 
of f ~ 0.3 a.u., which is remarkably independent of the angle of incidence in Fig. Ill While such a 
threshold could not be anticipated from the behavior of the average number of the lowest-energy vr 
plasmon mode j = 2, m = in Fig.lH it is perhaps noteworthy that this threshold speed correlates 
with the slope of the dispersion curve for the mode j = 2, m = when A; > 0.1 A^^ in Fig.lH 

In the limit of grazing incidence, ^ = 0, the particle travels parallel to the nanotube at a constant 
radial distance rmin. In that case, we can consider both internal (rmin < R) and external {r-^i^ > R) 
particle trajectories and define the stopping force acting on the particle, as described in the last 
paragraph of appendix |Bl Figure [8] shows the results for stopping force at different positions (a) 
inside and (b) outside a SWCNT. We observe that this force exhibits a peak structure which is quite 
similar to that observed in Fig.|7]for the total energy loss, with similar trends as the particle moves 
away from the nanotube wall. This can be rationalized by the fact that the shapes of the curves 



21 



for energy loss in Fig. |7] are largely independent of the angle of incidence, and by considering the 
stopping force as the energy loss of a particle per unit path length in the limiting case of oblique 
incidence when 6' — )■ 0°. 



IV. CONCLUSIONS 



The formulation developed in this work represents a direct but more realistic continuation of 



the works presented previously by us ||35i |46|] related to plasmon excitation by external charges 



moving paraxially in hollow cylindrical nano-structures. We applied a quantization procedure ll46ll 



to the two-fluid hydrodynamic model developed by Mowbray et al. 113211 . in order to obtain the 
average number of plasmons excited, and the total energy loss suffered by a fast charged particle 
passing near the surface of a single-walled carbon nanotube with an arbitrary angle of incidence. 

One of the most important results of this version of the hydrodynamic model is due to the 
inclusion of the restoring interaction, which causes the a+ix electron collective mode to oscillate at 
finite frequency in the limit of vanishing wavenumber, and changes the quasi-acoustic dispersion of 
the TT plasmon so that it agrees well with the available EELS data. On the other hand, quantization 
of these modes allowed us to obtain several quantities in terms of the average number of excited 
plasmons, such as the stopping force, energy loss spectra and total energy loss. We studied these 
quantities as functions of various parameters: the total speed of the incident particle, its minimum 
distance to the nanotube surface, and the inclination angle of its trajectory with respect to the 
nanotube axis. 

In addition, we have discussed (in appendix |B]) a relation between the quantized and semi- 
classical approaches. The latter approach presents the possibility of including the effects of plas- 
mon damping and hence enables a more direct comparison with the available experimental energy 
loss data, which is left for future work. 
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Appendix A: Statistics of plasmon excitations on carbon nanotubes 



In this appendix we use the methods of classical statistics to discuss energy losses of a fast 
charged particle due to electron plasma excitations on a carbon nanotube. We have seen that the 
quantization approach in the main text represents plasma excitations as a pool of independent 
quantum oscillators, or plasmons of mode q = (j, m, k) with an eigenfrequency Wq = cujmk- 



According to Carruthers and Nieto [|47ll . the probability of exciting N plasmons of any given 



mode q is given by a classical Poissonian distribution, 

Pq(iV) = e-'^'^ ^, (Al) 

where ;Uq = i^jmk is the mean number of plasmons excited in mode q. For further consideration, 
it is convenient to define the characteristic (or moment-generating) function associated with the 
probability Pq(iV), 

oo 

Gq(C) = (e^^^)q = e'^^ ^q(^) = e'^<'(^''-^). (A2) 

N=0 

It is sometimes of practical or theoretical interest to discuss marginal probability distribution 
Pjm{N) for having N plasmons excited in modes with given j and m values, while the longitu- 
dinal wavenumber k takes a full range of allowed values. Because of the statistical independence 
of plasmons of mode q, the excitation of plasmons with fixed (j, m) can be imagined as a sub- 
ensemble consisting of independent oscillators with different values of k. Then, the characteristic 
function associated with the marginal probability Pjm{N) can be written as a product of the char- 
acteristic functions G'q(C) for each member of this sub-ensemble, 

^,-(0 = \{G^{C) = n^'^''^*'''"'^ = e'^-(«'^-i), (A3) 

k k 

where = Xl/t /^q- ^y expanding the final result in Eq. (I A3 1) in a series of powers of the factor 
e*'', we find that the probability of exciting N plasmons in the mode (j, m) is also a Poissonian 
distribution, 

P,m{N) = e-^^- (A4) 

with the average number of plasmons, fijm, given in Eq. (|4T|) . 

Further, assuming that A'q is the number of plasmons excited in the mode q with frequency Wq, 
one can express the probability density P{£) for losing energy e as an ensemble average taking 
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into account the above mentioned statistical independence of plasmons of mode q, 

Pie) -- 



dr 

dr 
2jxh 



-ier/h 



exp 



q 



(A5) 



where the last step was derived by the use of Eq. (|A2|) . By expanding the final result in Eq. (IA5|) 
in a power series of the factor e*'^''^, we obtain the probability density as 

P(£) = 5^^q5(£-fe^q). (A6) 

q 

This can be written in a more explicit form by invoking Eq. (|36l) and using the delta function in 
Eq. (IA6I) to replace Wq = uojmk e/fi,a& 



P{e) — 6{e — hujmk) Cjmk 

jmk 



Kxt ( R; m, k, - 



E — / dk Kxt yR; m, ^ CjmkS {e - fhjjmk) 

n=-oo j=l 



(A7) 



where Cjmk is given in Eq. (l37l) . whereas Ven{R; m, k, cu) is given in Eqs. (l38l) . (l39l) . and (l40l) for 
a straight line trajectory. 

Finally, we note that the above result for probability density of energy loss can be used to 
evaluate the total energy loss as the mean value, A^Eioss = Jde e -P(e), so that 



2 oo „oo 
A-Eloss = / ^f^cifJ'Ci —^/^ / ^ — / dkUjmkl^jmk- 



(A8) 

ZTT / 

q j=l m=— oo 

It will be shown in the appendix |B] that the above results for P{e) and AE'ioss have a close relation 
with a semi-classical approach based on the hydrodynamic model. 



Appendix B: Semi-classical model of plasmon excitations on carbon nanotube 



A semi-classical, two-fluid 2D hydrodynamic model of plasmon excitations on a carbon nan- 
otube can be obtained by using the equation of motion Eq. ([3]) for the function ^y{vs^t) = 
z, t) augmented by the friction. 



m,. 



(Bl) 
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where 77^ > are friction coefficients for the z/th fluid, which can be used to describe plasmon 
damping in a phenomenological manner. 

The system of equations Eq. (IB II) . coupled with the Poisson equation Eq. © for the induced 
potential on the nanotube surface, can be solved by using the Fourier transform with respect to 
coordinates and time, {(p, z, t} — t- {m, k, to}, so that 



00 

iuit t 



iy{m,k,uj) = \ dte"^^^^{m,k,t), (B2) 
J —00 

with ^i/(m, k, t) defined in Eq. Q. Using the fact that the total induced number density of electrons 
per unit area is 5n = n° V/^o^ + r;,° V/.^^, we can write its Fourier transform as 

5n{m, k, to) = —x{m, k, u)Vext{R] rn, k, to), (B3) 

where the density response function of a carbon nanotube is given by 

y(m k u) = Xo{m,k,u;) 

with xo = + ' where the non-interacting response function of the z/th fluid is given by 

Sl [k^ + ^)+ul^-u{u + IT]^) 

It may be worthwhile quoting the final expression for k, to) in terms of the quantities defined 
in Eqs. (fTSi) . (fT6l) . and (|22|) and finite friction coefficients. 



X(m, k, u) = — — . ^ — t;^ — "—^ (k (B6) 

We define the total energy loss as the work done by the induced force on the external charge as 
it moves along its entire trajectory ro(t) with velocity vo(t) = ro(^), 

/oo 
dtvo(t) ■FiM(ro(t),t), (B7) 
-00 

where 

Fi„d(ro(t),t) = -eZ V$i„d(r, t)|r=ro(i) • (B8) 

On substituting Eq. (IB8I) in Eq. (IB7I) and using the chain rule, ^ = vq (t) • V + ^, we can eliminate 
the conservative part of the time integral in Eq. (IB7I) . thus arriving at 



Ai?ioss = -eZ I dt 



J —00 

00 7, . ^00 



J r=ro(t) 



^ u dte-'^'^indiro{t),u), (B9) 
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where $ind(r, co) is the Fourier transform of the induced potential with respect to time only, 

/oo 
da;e^-*<l>.„d(r,t), (BIO) 
-oo 

which can be obtained at arbitrary points outside the nanotube from 

$md(r,t) = — J ^e'"'''^"''9mk{r,R)5n{m,k,t), (Bll) 



m=—oo 



with gmk{r,R) = AnIm{\k\R)Kmi\k\r). 

Finally, by combining Eqs. (IB3I) . (IB9I) , (IBIOI) . and (IBl II) . and referring to the definitions in 
Eqs. (fT2)) and (|38l) . we can write the total energy loss as 



I /-oo ij^f^ ^ 2 

AE'ioss = ^ V / TTT-T? / c^t^t^ Kxt(i?;"i, /c,a;) 3 A;, w)] , (B12) 



where we have used the property that the real and imaginary parts of the density response, 
x(m, A;,u;), are an even and an odd function of frequency cu, respectively. Since the total en- 
ergy loss can be written as a first moment of a semi-classical version of the probability density for 
energy loss e = tku. 



A^ioss= / deePUe), (B13) 
Jo 

by comparison with Eq. (|B12|) . one can deduce 



(B14) 



It can be shown that the expression in Eq. (|B14|) reduces to the limiting expression in Eq. (IA7I) . 
which is obtained using the quantization of plasma excitations, i.e., when the plasmon damping is 
neglected. Namely, when both friction coefficients in the semi-classical hydrodynamic model are 
vanishingly small, rj,, — )■ 0+, we obtain for the imaginary part of the nanotube's response function 
at positive frequencies, u; > 0, 

3 [x(m, k, u)] = II^M J2 D^k S - ^jmk) , (B15) 

with the coefficient Dj^ak defined in Eqs. (I24l) . (|25] ) and (|26l ). Using Eq. (|37l) . this can be written 
as k, uj)\ = nhA J2j Cjmk 6 {uj — ^jmk), which, when substituted in Eq. (IB14|) . gives the 

probability density of energy loss Eq. (IA7I) of the quantization approach. Of course, unlike the 
quantization approach, it is possible, and even desirable to study the effects of plasmon damping 
on energy-loss distributions by using finite values for r/^ in the present semi-classical approach. 
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Moreover, the result in Eq. (IB 151) can also be used to show that the semi-classical result for 
the total energy loss, Eq. (IB 121) . is equivalent to the result Eq. (IA8I) of the plasma quantization 
approach, in the limit of vanishing plasmon damping. Of particular interest here is the case when 
the external perturbing charge moves parallel to a nanotube with speed v^^ (while v± = 0) at fixed 
distance tq > R, so that (po = and ZQ^t) = v^^t. We initially assume that the nanotube has a finite 
length L, and hence the traversal time of the external particle is T = L/v\\. Using the expressions 
Eqs. (fT2l) and (1381) . we obtain in the limit L — )■ oo, 

2 



Kxt (-R; m, k, u) = [e ZRgmk{R, tq 



T/2 

-T/2 
2 



^^g^(a;-fc^,||)^ 



27r— [e^ZRgmk{R, ro)] 5 [u - kv\\) 



(B16) 



allowing us to derive from Eqs. (IB 121) the average stopping force for the external particle, S 

AEioss/L, as 



POO 

S = m{e^Zfy^ / dkkll{\k\R)Kl{\k\r,)'^[x{mXkv\{)\ 



(B17) 



when ro > R, as shown in Fig.[8];b). On the other hand, when the perturbing charge is inside the 
SWNT, To < R, we obtain an analogous expression for the stopping force of 



S = m{e^Zfy^ / dkkll{\k\ro)Kl{\k\R)'^[x{mXkv\{)\ 



(B18) 



as shown in Fig. [8ta). It should be noted that these expressions for the stopping force are identical 



to those previously derived from semi-classical models of carbon nanotubes |132|, 
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Appendix C: Planar case 



We outline here how the formalism of the main text can be adapted to describe plasmon excita- 
tion in a planar 2D electron gas by an external charge moving on a specularly reflected trajectory 
based on the two-fluid model. There may be some interest for this outline in view of possible 
applications to EELS experiments on free-standing graphene M>M\- 

The planar case is simply retrieved by assuming that the nanotube radius grows indefinitely, 
i? — J- oo, in such a manner that the position along the nanotube circumference can be defined by a 
Cartesian coordinate y = R(p, whereas the radial distance from the nanotube wall can be defined 
by the Cartesian coordinate x = r — R. Similarly, with the longitudinal wavenumber k renamed 
kz, the ratio m/R becomes a quasi-continuous variable, ky, corresponding to a wavenumber for 
collective modes propagating around the nanotube's circumference. Consequently, the d(f integral 
and the m summation in Eqs. (|7]) and ^ become, respectively, 

TT CXD 

RjdLp---^ J dy-- - , (CI) 

— TT — OO 

OO ^ 

2nR ^ J 2n ' ^ ^ 

m=—oo 

— CXD 

indicating that we are using a 2D Fourier transform in the plane of the electron gas that maps 
{y,z} —7- {ky,kz} = K. Accordingly, the generic expression m? / R^ + k"^ appearing in the 
equations of the main text is mapped to K"^ = + ^^'^ the Green function becomes 

i?(7^fc(/2,r)^— e-^l^l, (C3) 

to the leading order in 1/ R. 

Next, we define a specularly reflected trajectory for the incident particle in Cartesian coordi- 
nates by ro(t) = {xo(t),?/o(^),-2o(i)} = {roit) - R,R(po{t),zo{t)} {x^m + V±\t\,Viiyt,Vii^t}, 
where Xmm > is the minimum distance from the planar electron gas and V_i_ is the normal com- 
ponent of the projectile velocity. Then, Eq. (fT2l) gives in the limit R oo 

Kxt(i?; m, k, t) ^ -e^Z f^e-^(--'"+^^l*l)-*K-V|i*, (C4) 

K 

where Vy = {V\\y, V\\z} is the parallel component of the projectile velocity. As a consequence, Eq. 
(1381) gives 

tAR;m,k,co) ^ -e^Z^e-""^--- ^ -, (C5) 
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so that the average number of plasmons in the mode (j, m, k) (j, K) finally follows from Eq. 
(I36l) as 



2no /^|k 



(a;,K-K-V||)' + (KV^j_ 



2 ' 



(C6) 



where Dj^ is given by Eqs. (|24l) . (|25] ) and (|26l ). in which plasmon frequencies are to be used from 
Eqs. (fT5l) . (fT6l) and (|20|) upon the replacements 

(C7) 



ml 



and 



^0 ^0 



(C8) 



It is worth mentioning that, if one lets both u^r — )■ and s,^ — )■ in Eq. (IC7I) . the planar 
version of Eq. (l20l) no longer gives rise to plasmon splitting, but rather one recovers (jJ2k — ^ and 
i^iK — ^ ^Txe^—K describing the familiar plasmon dispersion of a single-fluid model for planar 



2D electron gas with surface density uq = n '^ + and an effective mass [!17!]. 
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